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^ , Abstract 

2 ■ We consider the galactic population of gamma-ray pulsars as possible sources of 

cosmic rays at and just above the "knee" in the observed cosmic ray spectrum at 
IQiS-iQis eV. We suggest that iron nuclei may be accelerated in the outer gaps of 
• pulsars, and then suffer partial photo-disintegration in the non-thermal radiation 

fields of the outer gaps. As a result, protons, neutrons, and surviving heavier nuclei 
O I are injected into the expanding supernova remnant. We compute the spectra of 

^ ' nuclei escaping from supernova remnants into the interstellar medium, taking into 

■ account the observed population of radio pulsars. 

Our calculations, which include a realistic model for acceleration and propaga- 
tion of nuclei in pulsar magnetospheres and supernova remnants, predict that heavy 
nuclei accelerated directly by gamma-ray pulsars could contribute about 20% of the 
observed cosmic rays in the knee region. Such a contribution of heavy nuclei to 
the cosmic ray spectrum at the knee can significantly increase the average value of 
(In^) with increasing energy as is suggested by recent observations. 



PACS Numbers: 98.70. Sa (Cosmic rays), 97.60. Gb (Pulsars). 

Keywords: cosmic rays, spectrum and mass composition, pulsars, cosmic ray acceleration 

1 Introduction 

The cosmic ray spectrum is usually described by a superposition of power laws with in- 
dex ~ 2.7 below and ~ 3.2 above the so-called "knee" at a few lO^^eV. In the region 
of 10^^ — 10^^ eV the spectrum flattens, and has an index again close to 2.7. There are 
also suggestions that the shape of the spectrum in the knee region has a fine structure 
which may reflect acceleration of specific groups of nuclei, e.g helium, oxygen, and iron. 
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prominent in the Galactic abundances (see [|T], and references therein). Recent exper- 
imental results show that the cosmic ray composition changes through the knee region 
of the spectrum, becoming heavier with increasing energy [0, ^, 0. For example, the 
CASA-MIA experiment measures the mass composition in the energy range 10^^ — 10^^ 
eV [0, and shows that above 10^^ eV the iron group nuclei should dominate. The analysis 
of hybrid data from the High Resolution Fly's Eye prototype and the MIA muon array 
indicates that the mass composition becomes lighter again in the region 10^'' — 10^^ eV 
0, and is proton dominated at around 10^^ eV [0]. 

It is usually argued that cosmic rays with energies below 10^^ eV are accelerated 
at shocks of expanding supernovae (see e.g. ref. It seems that the acceleration 

mechanism of the higher energy cosmic rays may be different, although supernova shocks 
of various types have also been proposed B and also lead to fine structure in the spectrum 



and composition in the knee region [|T0|. This idea is consistent with the recently observed 
variation in the composition of cosmic rays in the region 10^^ — lO^^eV mentioned above. 
However, it is uncertain what fraction of the cosmic rays above 10^"^ eV, are actually 
accelerated by supernova shocks, and so other mechanisms able to accelerate cosmic rays, 
particularly heavy nuclei, to energies above lO^^eV may be required. It has also been 
suggested that protons and heavy nuclei could be accelerated to such energies in pulsar 



winds 13, 14 



Nuclei can be also accelerated in pulsar magnetospheres as discussed by Cheng, Ho 
& Ruderman ||T5|, Cheng & Chi |T6[ and Bednarek & Protheroe [|17|. In our previous 
work ||T7| , |l^, we considered acceleration of iron nuclei in the magnetospheres of very 
young pulsars, and Crab-type pulsars, in the context of 7-ray and neutrino emission by 
supernova remnants. We found that, although the photo-disintegration process of iron 
nuclei during their acceleration and propagation through the neutron star magnetosphere 
is important, heavy nuclei would be injected into the medium of expanding supernova 
remnant. In the present paper we compute the possible contribution of nuclei injected by 
the galactic population of gamma-ray pulsars to the observed cosmic ray spectrum. We 
show that gamma-ray pulsars can contribute to the cosmic ray spectrum with energies 
above ~ 10^^ eV. 



2 The pulsar spin evolution 

Neutron stars are likely to be born with initial periods Pq < 10 ms if we assume that the 
angular momentum of the pre-supernova star is conserved during the collapse. However, 
the initial angular momentum may be lost in the case of hot rapidly rotating neutron 
stars by emission of gravitational radiation due to the r-mode instability [|19| , 0, ^ 



As a result, the star would be left with a period of about 10 ms at about 1 year after its 
formation. After that, the period of the neutron star would change, mainly as a result of 
emission of magnetic dipole radiation with power given by 

= B^Rf^s^'^ sin^ i/6c^ (1) 
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where B is the surface magnetic field at the magnetic poles, i?NS is the radius of the 
neutron star, i is the angle between rotational and magnetic axes, and Q = 27t/P, where 
P is the period of the star. Also, gravitational radiation may be important for the case 
of stars having significant ellipticity k, 

Lgr = 32G'/^A;2fiV5c^ (2) 



where G is the gravitational constant, / is the moment of inertia |]23 . 

We have investigated the population of young pulsars observed with age < 10^ years 
(about 100 pulsars) from the catalogue of Taylor, Manchester & Lyne [2l| in order to 



estimate the neutron star ellipticity. By starting with an assumed initial period Pq = 
10 ms, perpendicular component of magnetic field Bsini, and ellipticity k we may use 
Eqs. |l] and ^ giving the rate of loss of rotational kinetic by magnetic dipole radiation and 
gravitational radiation, respectively, to calculate the period P and period derivative P as 
a function of time since birth. For each of the young pulsars considered, and assuming 
standard values of / and -Rns; and the inferred value of Bsini, we obtain the value of 
k such that P and P match those observed p5[. We note that the values of k obtained 



are rather insensitive to the assumed initial period for pulsars older than ~ 10^ years. 
The results are plotted in Fig. 1, and show that a pulsar's ellipticity must be correlated 
with the perpendicular component of its surface magnetic field multiplied by its radius 
and divided by its moment of inertia. Since the theoretical radii and masses of neutron 
stars lie in a rather narrow range, we conclude that there exists a correlation between 
the ellipticity and the surface magnetic field for neutron stars. The correlation coefficient 
obtained is equal to ~ 0.65, which for 100 pulsars means that the correlation is highly 
significant. The inferred dependence of pulsar ellipticity on surface magnetic field can be 
well approximated by 

log k = -22.3 + 1.571og B sin i, (3) 

where B is in Gauss. 

The above result suggests that the ellipticity of a pulsar may be due to the deformation 
of the shape of the neutron star by its strong interior magnetic field (see e.g. [^). The 



prescription for pulsar spin evolution we apply in our present work is based on Eqs. ^ 
Note, however, that in the case of some specific pulsars, the difference between the result 
from Eq. 3 and estimates based on the known age of the pulsar is significant, e.g., in the 
case of the Crab pulsar for which the estimated ellipticity is ~ 3 x 10~^ [^. In the case 



of the recently-discovered pulsar in the historical (AD 386) SNR Gil. 2-03, the standard 
formula |2^ gives for the pulsar B = 1.7 x 10^^ G, and using the historical age t = 1612 
yr a pulsar initial period of 62 ms is found |^ . If gravitational radiation were included, 
the initial period could be much shorter. We find that for the same magnetic field, an 
ellipticity k = 1.49 x 10~^ would give an initial period of 10 ms for AD 386. This value of 
k, together with B = 1.7 x 10^^ G, is consistent with the correlation described by Eq. 3 
and seen in Fig. 1. 
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3 Acceleration of nuclei by young pulsars 



Periodic emission of photons from radio to 7-ray energies gives direct evidence of the 
existence of electric potentials in neutron star magnetospheres which are able to accelerate 
particles to energies of at least ~ 100 GeV. To explain these observations, a number of 
models have been proposed including the polar gap model, the slot gap model, and the 
outer gap model (see e.g. Arons [Q). In the present paper, we concentrate on the 
outer gap model which is more attractive from the point of view of acceleration of nuclei 
by pulsars because it postulates induction of larger electric potentials in the pulsar's 
magnetosphere. The inner gap (polar cap) models (e.g. Ruderman & Sutherland 
Harding & Muslinov |^) predict that the electric potentials induced in the inner gaps 
are significantly lower (of the order of ~ 10^^ V), than expected in the outer gaps for 
pulsars with the same parameters. Nuclei accelerated through such potentials would not 
suffer efficient photo-disintegration in collision with thermal photons with a characteristic 
temperature ~ 10^ K emitted by the polar cap or from the whole surface of a classical 
radio pulsar (only very young pulsars with ages of days to a month are sufficiently hot to 
photo-disintegrate Fe nuclei [@). Hence, we neglect this stage of pre-acceleration of Fe 
nuclei during their passage through the inner gap. Moreover, the inner and outer gaps do 
not need to extend along these same magnetic field lines, and so they do not need to be 
physically connected. 

Nuclei, most probably iron nuclei, can be stripped off the hot surface of a neutron star 
due to the strong electric fields and the bombardment by particles produced in the pulsar 



magnetosphere [^, ^ . This can occur where the magnetic and rotation axes are counter- 
aligned, and therefore where the co-rotation charge is positive above the poles, and can 
pull Fe nuclei from the surface. If one assumes that pulsars have random orientations of 
these axes, then Fe nuclei can be extracted from a significant fraction of pulsars, but if 
this is not the case then our results could be significantly affected. 

These nuclei are then accelerated by the pulsar outer gap electric field to high energies. 
We assume that the process of extraction of iron nuclei from the pulsar surface occurs 
during the 7-ray emitting stage of the pulsar as the production of 7-rays in the outer gap 
requires efficient production of relativistic plasma, and this also heats the polar cap 
region. Chen & Ruderman |^ distinguish two types of gamma-ray pulsars, so called 
Crab-hke (Crab- type) and Vela-like (Vela-type). Beyond the Vela-like death line the 
gap cannot close, and this prevents copious production of secondary pairs, and causes 
efficient production of gamma-ray photons to cease. However, electric potentials inside the 
outer gap remain, and the acceleration of primary particles is still possible. In more recent 
papers (Romani Zhang & Cheng Cheng & Zhang [^) a third type of pulsar, so 
called Geminga-like pulsars, is distinguished. In these pulsars the production of secondary 

pairs inside the outer gap is sustained by the soft X-ray emission from the polar cap 
region. The death line for Geminga-like pulsars occurs when there is insufficient fiux of 
thermal photons from the polar cap region of the pulsar. However, primary particles can 
be still accelerated inside the outer gap even after a pulsar has crossed the Geminga-like 
death line (only pair production stops). We expect that the injection of the Fe nuclei 
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from the region of the polar cap is related to the thermal radiation field produced by this 
region of the neutron star surface. The thermal radiation field drops significantly at a 
neutron star age ~ 10^ yr (see Fig. 2 in Romani ||3^). We estimate that at this age the 
pulsar stops injecting Fe nuclei into the magnetosphere. 



3.1 Model for the electric field in the outer gap 

During the 7-ray production stage of pulsars, three different types of outer gap can be 
distinguished, depending on the relative importance of different processes of 7-ray produc- 
tion: the Crab-type and Vela-type gaps (Cheng, Ho, Ruderman ll5|), and the Geminga- 
type gap (Romani [^, Zhang & Cheng The separation between the Crab-type and 

the Vela-type outer gaps depends on the pulsar period and surface magnetic field, and is 
given by Eq. 25 of Chen & Ruderman ||3^ 
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(4) 



where B12 is the pulsar surface magnetic field in 10^^ G, and the period Pcv is in seconds. 
The maximum potential difference across different types of outer gap can be approximated 
by the simple formula (Cheng ||36 



private communication) , 



^{B,P) ^ 6.6 X 10'' f'BuP-' V. 



(5) 



This formula gives the potential drop in outer gaps a factor of two lower than derived in 
original paper by Cheng, Ho & Ruderman T5[| , 



and is more consistent with the recent 



works on the outer gap model [^, ^ 



However, we note that this simple formula may not 
match all observational data well because of uncertainties of viewing angle and inclination 
angle which are very important in determining the size of the outer gap. The parameter 
/ is used to parametrize the gap voltage and power, and is defined by Cheng, Ho & 
Ruderman IITS 



/ ^ 825^2''/'°^''/'°, (6) 

Formulae (5) and (6) give the potential difference ~ 3 x 10^^ eV for the Crab pulsar. For 
Geminga-type gaps (Zhang & Cheng [ p4[ ) 



/ ^ 5.5Sr2'/V26/2i. (7) 

For Vela-type outer gaps a simple definition of the parameter / is not available. We 
assume that in this case the potential drop in the outer gap is equal to a constant fraction 
of the total potential drop in the outer gap given by Eq. (5) with p = 0.15 for all periods 
^). The dependence of the potential difference across the outer g function 



see 



of the pulsar period is plotted in Fig. g for three different values of the surface magnetic 
field of the pulsar. 
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According to Cheng, Ho & Ruderman O] the potential drop across the gap depends 



also on the position of injection into the outer gap field, z (z ranges between and a), in 
relation to the maximum gap width a. Therefore, we assume that the potential drop for 



position z is given by (see definitions below eq. 3.10 in ref. ||T5| ) 



$(z) =4<I>(5,P)^[1-^]. (8) 
3.2 Radiation field inside the outer gap 

Iron nuclei accelerated in the outer gap potential can suffer significant photo-disintegration 
in collisions with photons produced inside the gap. We estimate the density of soft photons 
inside outer gaps of pulsars with a given period and magnetic field by scaling from the 
density of soft photons inside the outer gap of the Crab pulsar assuming the outer gap 
model. This density can be estimated from the observed pulsed fiux of soft photons, 
assuming that the dimension of the outer gap is of the order of the light cylinder radius. 
Ho 1^ has shown, by analyzing a sample of Crab-type pulsars with periods ~ 15 ms, 
that the luminosity of the photons produced in the gap does not change drastically for 



Crab-type outer gaps. Based on the results of Ho [^, we interpolate/extrapolate the 
density of photons in the present Crab outer gap for other periods. The 7-ray luminosity 
of pulsars scales as ~ times the total spin down pulsar luminosity [|1^, and so the 
7-ray luminosity of the Crab type pulsars depends on the magnetic field as ~ iJ^/^o (gee 
Eq. 1^), i.e. it is only very weakly dependent on B. Since energies of synchrotron photons 
are proportional to the magnetic field in the gap, then the photon density in Crab-type 
outer gaps should be inversely proportional to the pulsars' magnetic fields. 

We note that the observations seem to suggest that the 7-ray luminosities of different 
types of pulsar are proportional to BP~^ (e.g. Fig. 7 in Arons [^]), although different 
models predict different dependencies of 7-ray luminosities on B and P [Q. However the 
above dependence should be considered with some caution because it interpolates between 
the Geminga pulsar (and PSR 1055-52) and the Vela pulsar (and its twin PSR 1706-44). 
The pulsar PSR 1951-1-32 deviates by an order of magnitude and the Crab pulsar by a 



factor of two, and if the distance to the Vela pulsar is ~ 250 pc but not ~ 500 pc 
then Vela pulsar deviates by a factor of four from this dependence. Also, the distance to 
Geminga pulsar is not well known. 

The density of soft photons in Vela-type outer gaps is obtained based on the observed 



change of the pulsed luminosity of the Vela type pulsars (e.g. [^), and our estimate of 
/ in the case of Vela-tj^e pulsars, with the normalization to photon density in the Crab- 
type outer gap at the transition period given by Eq. (^. A good fit to the luminosity of 
the Vela-type pulsars is obtained if the density of photons in the outer gap drops with 
period as oc P~^. Therefore, we apply the following model for the soft photon density in 
Crab and Vela-type pulsars for different periods and surface magnetic fields 

(9) 




Pora^gorah ^ tj^g Crab-type; 

Pc.abgc.ab (^P(^^ ^ ^ ^i^g Vela-type, 
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where Pcrab = 0.033 seconds, -Bcrab = 4 x 10^^ G and is the density of photons in 

the outer gap of the Crab pulsar at present (Cheng, Ho & Ruderman 1986). 

The Geminga-type outer gap is determined by the thermal photons produced on the 
neutron star surface. This thermal radiation field for pulsars with ages corresponding to 
Geminga-type pulsars is too low to create a significant target for photo-disintegration of 
Fe nuclei. The epoch of gamma-ray emission of Geminga-type pulsars is determined by 
the temperature of the neutron star surface which, according to models of neutron star 
evolution, drops rapidly after ~ 10^ years (see e.g. Fig. 2 in Romani Therefore, 
Geminga-type outer gaps should not be able to operate after ~ 10^ years from the ex- 
plosion of supernova. Here we assume that the process of extraction of Fe nuclei from 
the neutron star surface does not occur after this age because, as we noted above, it is 
probably caused by heating of the polar cap region by pairs created in the pulsar 
magnetosphere. 



3.3 Number of nuclei accelerated in the outer gap 



As in our previous paper [17], we assume that the number of Fe nuclei injected per unit 



time, A^Fe; is simply related to the total power output of the pulsar Lem(-B, P) IH3i , 

iVpe = e^em(5, P)/Ze<^{B, P) , (10) 

where ^ is a free parameter describing the fraction of the total power taken by relativistic 
nuclei accelerated in the outer gap, Z = 26 is the atomic number of Fe, and $(-B, P) 
is defined by Eq. |^. The parameter ^ has to be always less than one since the power 
of accelerated Fe nuclei can never be greater than the spin-down power of the pulsar. 
The charge density estimated from Eq. [1^ exceeds the 'classical' Goldreich- Julian charge 
density close to the polar caps p = —{QB /27ic)[l — {fir / cY sin"^ 9]~^ ~ —QB/2ttc, where 
the approximate formula is valid close to the surface where Qr sin 6 <^ c. However, as 



noted by Goldreich & Julian below their Eq. (8), the above formula is only valid 
in the co-rotation portion of the magnetosphere. If there is additional negative current 
onto the surface accompanying the charge injected from the surface in the form of Fe 
nuclei, then the above formula is not valid. Positive charge may flow onto the surface 
along the same magnetic lines as Fe nuclei, and could be produced somewhere in the 
pulsar magnetosphere, not necessary close to the pulsar's surface (e.g. possibly in the 
outer gap). Pair production just above the polar cap does not need to occur in this case. 



3.4 Number of neutron stars with different parameters in the 
Galaxy 

As discussed above, we assume that all pulsars are born with periods of about 10 ms, or 
that they reach this period soon after formation by emission of gravitational waves due 
to the r-mode instability. As we argued in Section 2, pulsars need to have a significant 
ellipticity k in order to reach their observed periods from a period at birth of 10 ms. 
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The distribution of the pulsar magnetic fields can be evaluated from observations, 
assuming that the surface magnetic field does not decay significantly with pulsar age. For 
example, Narayan [J5| approximates the distribution of surface magnetic fields of pulsars 



at birth by 

dN/d{\ogB) ^ 0M5{W^^G/B) yr^\ (11) 

for i? > 2 X 10^^ G. We estimate the distribution of pulsars with surface magnetic fields 
below 2 X 10^^ G based on the Fig. (13) in Narayan [^. It is likely that this approxi- 



mation underestimates the number of pulsars with relatively low magnetic fields because 
of observational selection effects. The above formula is normalized to a formation rate of 
pulsars in the Galaxy equal to 1 per 70 years. 



4 Acceleration and propagation of nuclei from the 
neutron star surface 

Nuclei extracted from the neutron star surface are accelerated and photo-disintegrated 
in the outer gap. Their disintegration products, and those which escape without any 
disintegration, are injected into the nebula and lose energy in collisions with matter ejected 
during the supernova explosion, and as a result of the expansion of nebula. All these 
processes are discussed below. 



4.1 Propagation of nuclei through pulsar magnet ospheres 

The importance of photo-disintegration of nuclei with different mass numbers during 
propagation of nuclei in the radiation field of the outer gap is determined by the reciprocal 
mean free path, which we calculate using cross sections of Karakula & Tkaczyk |^ (with 
improvements given in ref. [|T3]). The results of such a calculation for an outer gap with 



parameters applicable to the Crab pulsar are shown in Fig. 1 of Bednarek & Protheroe [|17 
It is evident that multiple photo-disintegrations of primary Fe nuclei will occur for this 
pulsar. In our present calculations, the photo-disintegration of nuclei in the thermal 
radiation field emitted by the pulsar surface is neglected since it may be only important 
during a few days after creation of the pulsar when its surface temperature is of the 
order of 10^ K and consider only photodisintegration in the outer gap's non-thermal 
radiation field. 

We simulate the photo-disintegration process of nuclei accelerated by pulsars using 
the evolution model, initial parameters of pulsars, the outer gap model for acceleration, 
and the radiation field discussed in Sections 2 and 3. Multiple photo-disintegrations occur 
only for pulsars with Crab and Vela-type outer gaps. We simulate the acceleration and 
propagation of iron nuclei in pulsar outer gaps in order to obtain the energy spectrum of 
neutrons and protons (extracted from Fe nuclei) and the spectra of escaping secondary 
nuclei. Note, however, that the fate of neutral particles (neutrons) and charged particles 
(protons and secondary nuclei) is different. The neutrons move ballistically through the 
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surroundings of the pulsar, and decay at distances determined by their Lorentz factors, 
and this can be inside or outside of the expanding envelope of the supernova. We compute 
the spectra of neutrons (and so protons from their decay) which decay inside and outside 
the Nebula as in Bednarek & Protheroe . In the case of protons from neutrons decaying 
outside the nebula, we assume that all of them escape into the Galactic medium without 
being overtaken by the expanding nebula (we drop the second term in brackets in Eq. (5) in 
Bednarek & Protheroe Jl^ which represents the effects of protons from neutrons decaying 
outside the nebula being overtaken and trapped by the expanding nebula at some later 
time). This term can be neglected since the minimum diffusion distance of these protons, 
with Lorentz factor 7p = 10^, in the Galactic magnetic field 5 x 10~^ G is comparable to 



the radius of nebula during its time of free expansion [|T7 



Protons and secondary nuclei injected into the Nebula are trapped by the magnetic 
field of the nebula, and suffer adiabatic and collisional energy losses during the time of 
expansion of the nebula. These effects are considered in the next subsection. 



4.2 Interaction of nuclei with matter inside nebulae 

Soon after the supernova explosion, when the nebula was relatively small, nearly all the 
energetic neutrons would be expected to decay outside the nebula. However, at early 
times we must take account of collisions with dense matter of the expanding nebula. The 
optical depth for neutrons moving ballistically through the envelope may be estimated 
from ~ (ynpTLijr ^ 8.6 X lO^'^Mif^^t"^, where M = MiMq is the mass ejected during 
the supernova explosion in units of solar masses, r = vt, = M/(4/37rr^mp) is the 
number density of target nuclei, and v = lO^fg cm s~^ is the expansion velocity of the 
nebula. We note that TnH = 1 at t ^ 0.93Ml^^Vg^ yr. 

Charged nuclei (protons, protons from decay of neutrons, and secondary nuclei) in- 
jected by the pulsar are captured by the nebula's magnetic field and interact with the 
dense material of the expanding supernova remnant. The optical depth for these nuclei 
depends on the time, t, of the injection of nuclei, and can be estimated from 

rAp(t) ^ 1.3 X lO^'^A^/^Miv^h-^ (12) 

Interactions of nuclei with matter inside the nebula become negligible when the optical 
depth becomes less than one. For the parameters derived for the Crab, f § = 2 and 
Ml = 3, and the average value of the mass number of nuclei injected from the pulsar 
magnetospheres into the nebula < A >pa 48, the nebula becomes transparent for nuclei 
at times t > 3.6 x W^A^/^mI^^v^^^^ ^ 22 years. Nuclei injected at earlier times should 
suffer significant fragmentation in collisions with the matter inside the nebula. We have 
simulated this process of fragmentation of nuclei injected by the pulsar applying the semi- 
empirical formulae for the cross sections for the fragmentation of nuclei of Silberberg & 
Tsao with parameters given in their tables 1A,B,C,D. In our simulation, we assumed 
that in a single collision a nucleus fragments into two nuclei with different mass numbers. 

Protons created during complete fragmentation of nuclei are considered as being cap- 
tured by the magnetic field of the nebula. In contrast, neutrons which move ballistically 
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through the nebula have some chance of decaying inside or outside the nebula. These 
effects have been included in our simulation code. We consider the interactions of these 
secondary protons and neutrons with matter inside the nebula in the same way as for pro- 
tons and neutrons from photo-disintegration of iron nuclei in the pulsar's magnetosphere. 

4.3 Escape of nuclei from nebulae 

We assume that nebulae expand during the free expansion phase with the velocity ob- 
served in the case of the Crab Nebula, i.e. 2000 km s~^. When the amount of matter 
ejected during the supernova explosion becomes comparable to the amount of interstellar 
matter swept up by the nebula, then the nebula reaches the Sedov phase. The Sedov 
phase starts when the nebula reaches radius 

Rs = 6.8 X 10^^(Mi/n)^/^ cm, (13) 

and this occurs at time 

ts = Rs/v, (14) 

after the explosion. Here n is the density of the surrounding medium (cm~^), and v is 
the expansion velocity of supernova during its phase of free expansion. For a supernova 
with an ejecta mass of Msn = 3M0, a velocity of free expansion v = 2000 km s~^ as 
observed for the Crab Nebula, and a density of the interstellar medium n = 0.3 cm~^, 
the Sedov phase should begin when the nebula reaches a radius of i?s ~ 4.8 pc at about 
ts ~ 2000 years. This is only about twice as large as the present radius (and age) of the 
Crab Nebula. We assume the following approximation for the radius of a nebula at time 
t 



^Neb(t) = {j;|/5^2/5^ Ih^J ^^^^ 

The magnetic field inside the nebula changes with time and depends on the distance 
from the center of the nebula, i.e. the pulsar. We estimate it applying the results of the 
Kennel & Coroniti model, which gives the distribution of the magnetic field inside 



the nebula as a function of the parameter a (the ratio of the magnetic energy flux to 
the particle energy flux at the location of the pulsar wind shock at radius -Rsh)- We 
estimate the magnetic fleld inside the nebula (in units of the magnetic fleld strength Bgh 
at the shock location) as a function of its radius based on the flgs. 3 and 7 in Kennel & 
Coroniti |^8|]. For nebulae with cr ^ 0.003, the magnetic fleld drops signiflcantly with the 
distance from the center of the nebula (flg. 3 in ||4^), and can reach a value typical for 
the interstellar medium (5 x 10~^ G) at distance which is less than the radius of the 
expanding nebula -RNeb, and so we assume that beyond radius -Rb, the magnetic fleld is 
equal to 5 x 10"^ G. 
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The magnetic field at the shock, Bsh, can be found from the dependence of the pulsar's 
magnetic field in the region of the inner magnetosphere {B oc R"^), and the pulsar wind 
zone {B oc R~^). It is given by 




where -Rs and B^ are the radius and the surface magnetic field of the pulsar, and i?L = 
cP/2'K is the light cylinder radius of the pulsar. The value of a is estimated for the Crab 
pulsar to be ~ 0.003 (e.g. Kennel & Coroniti |4^). However, it has recently been found 
that for the Vela pulsar the value of a is rather close to 1 (Helfand et al. |^). To our 
knowledge, there is no available information on a for other plerions. Therefore, in order to 
fulfill observational constraints of these two pulsars, we assume that pulsars with periods 
lower than the period of the Crab pulsar have a ~ 0.003 and pulsars older than the Vela 
pulsar have a ~ 1. For pulsars having a rate of the rotational energy loss between that of 
the Crab and Vela pulsars, a is obtained by linear interpolation between these two values. 

The location of the wind shock, i?sh, can be estimated from the comparison of the 
wind energy flux at the shock with the pressure of the outer nebula, which is determined 
by the supply of the magnetic energy to the nebula by the pulsar over its all lifetime (Rees 
& Gunn [0). The location of the shock at an arbitrary time tp after explosion of the 
supernova can be estimated from 



-^cm(^p) ^ So' 0'Lem{t)dt , , 

iTiRlc ^ InRU ■ ^ ' 

As an example, we show in Fig. 3 the characteristic radii in the nebula, -RNeb, Rsh, 
and Rb, for the following parameters: Bg = 3 x 10^^ G, A; = 3 x 10""^, v = 2000 km s~\ 
^SN = 3M0, and n = 0.3 cm~'^. The location of the shock, estimated from the above 
formula, is consistent with the observations of the shock in the Crab [HSI and Vela [EU 



nebulae. For older nebulae Rb becomes less than -RNeb, i-e. the magnetic field inside the 
outer regions of expanding nebula drops to 5 x 10~^ G. 

Charged particles, nuclei and protons, injected into the nebula will diffuse, and finally 
escape from it into interstellar space. We assume that particles escape from a nebula at 
time tinax when the diffusion distance of particles injected at the time t becomes compa- 
rable to the radius of the nebula -RNeb at the time tmax- In order to estimate this time at 
which they escape from the nebula, we compare the diffusion distance for nuclei with the 
mass number A and the Lorentz factor 7a in the time dependent magnetic field inside 
an expanding nebula with radius given by Eq. The distance diffused depends on the 
diffusion coefficient which we take to be 

D = D^p"\ (18) 

as appropriate to a Kolmogorov spectrum of turbulence. Here, p ~ E/eZ is the particle 
rigidity, and Dq is normalized in such a way that D = ri^c/3 when the Larmor radius. 
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tl, is equal to the radius of the remnant (largest scale) at time t (given by Eq. |T5|). This 
is a conservative assumption as it will give rise to the greatest trapping of nuclei in the 
nebula (and highest adiabatic losses) while being consistent with a Kolmogorov spectrum 
of turbulence. However, the precise form adopted for D has little effect on our results as 
most particles are released at tmax- The diffusion distance is obtained by integration over 
time from the injection of nuclei at tinj up to escape at tesc 




^dif = / \hlidt'. (19) 



Note that the diffusion coefficient depends on the age of the nebula, and this is included 
when integrating Eq. |19|. For nuclei with sufficiently low energies, tmax may be shorter 
than tesci and we then take tesc = ^max- The adiabatic losses of particles inside nebula are 
considered up to the moment when the expansion velocity of the nebula during the Sedov 
phase, fexp = v{ts/ty^^, drops to a few tens km (i.e. to the characteristic velocities 
of the clouds in the interstellar medium). This happens at the age of ~ 10^ yr. Note that 
at this age the magnetic field inside the nebula is larger than 5 x 10^^ G only at small 
region above the radius of the shock {R^ ~ 2i?sh, see Fig. 3). 

Adiabatic and interaction energy losses of nuclei during the expansion of the nebula 
are included by estimating the Lorentz factor of nuclei at arbitrary time, t, 

where K is the inelasticity coefficient, assumed equal to 0.5 for protons and neutrons. For 
nuclei with mass number A, we take the inelasticity coefficient to be i^' ~ 0.5 /A^^^ since 
the probability of a collision of a single nucleon within the energetic nucleus, with a proton 
of the remnant nebula is, on average, proportional to (aAp/A) oc (A^^^/A) oc A~^^^. 



5 Injection rate of nuclei into the Galaxy 

Taking into account all of the effects discussed above, we compute the injection rate of 
nuclei of different mass numbers by the galactic population of pulsars into the interstel- 
lar medium, and this is given in Fig. ^(a). The spectrum shows two distinct maxima 
corresponding to heavy nuclei (dot-dot-dot-dash curve peaking at a few lO^^eV) and to 
protons (solid curve peaking at ~ lO^^eV). The contribution of medium mass nuclei is 
visible at energies between the proton maximum and the heavy nuclei maximum. The 
main contribution to the spectrum of nuclei below the high-energy peak is from Crab-type 
pulsars. This is clearly seen in Fig. |^(b) in which we present the contribution of pulsars 
having different types of outer gap to the total spectrum of nuclei injected into the Galaxy 
the from pulsars. The main contribution to the maximum at 3 x 10^ GeV comes from 
Vela-type and Geminga-type pulsars. However the nuclei with the highest energies are 
injected by the Vela and Crab-type pulsars as expected from inspection of the electric 
potentials shown in Fig. ^. 
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6 Cosmic ray propagation 



We assume that the radial distribution of the cosmic ray emission from gamma-ray pulsars 
follows the distribution of old neutron stars in the disk of the galaxy as given by Paczynski 



51 



q{r)2'Krdr = aRex.Y>{-r / Rcy^^)Rj^2'nr dr (21) 

with aji = 1.0683/27r, i?cxp = 4.5 kpc, and r being the the distance from the centre of 
the Galaxy of radius rmax = 20 kpc, and that they are uniformly distributed in a disk of 
thickness 2hg. For hg we take an estimate of the scale height of pulsar birth locations, 
hg = 130 pc since random pulsar velocities ~ 200 km s"^ would not significantly 



increase the scale height for pulsars much younger than 10^ y. Combining q{r) with the 
calculated spectra of cosmic rays of atomic number Z, Qz{E) = dN/dE/dt presented in 
Fig. §(a), we obtain the spatial distribution of cosmic ray injection 

ar-E) = qir)QziE)/{2hg). (22) 

We assume that at energies above 10'^ GeV we can reasonably neglect interactions 
and energy losses of cosmic ray nuclei, and that diffusion is the dominant process in 
shaping the observed cosmic ray spectrum. We consider two different models (Ptuskin 
et al [^) for propagation of cosmic rays in the galaxy which are consistent with the 
observed secondary to primary ratios: (i) the minimal reacceleration model, and (ii) the 
standard diffusion model. We use the best-fitting parameters found by Ptuskin et al 
for a halo half-thickness hh = 5 kpc (we note here that the cosmic ray intensity at Earth 
obtained is rather insensitive to the value of hg chosen provided that hg -C hh). For the 
minimal reacceleration model, the diffusion coefficient used is 



D = 3.8 X 10'^(p/l GV)^/^ cm's"^ (23) 

where p ~ E/Z is the particle rigidity. For the standard diffusion model, the diffusion 
coefficient used is 

D = 3.0 X 102^(p/l GV)°-^' cm^s-^ (24) 



Assuming cylindrical geometry, we then use Eq. 3.22 of Berezinsky et al. [Q, to calculate 



the cosmic ray density at Earth assumed to be in the Galactic plane 8 kpc from the centre 
of the Galaxy. 

The spectra after propagation are plotted in Fig. | for the various species where they 
are compared with the world data set of observed cosmic ray intensity [Q. Fig. |^(a) and 



(b) show results for the two propagation models. The results are rather insensitive to the 
propagation model used, and this can be easily understood after inspecting the injection 
spectra in Fig. H(a) and Eqs. [23|-^. The peak at 10^ GeV is due to protons with rigidity 



10^ GV, and the peak at ~ 3 x 10^ GeV is due to iron nuclei with rigidity 1.5 x 10^ GV, 
approximately the same as the proton peak. Looking at Eqs. ^ and |2^, we note that the 
diffusion coefficients for the two propagation models are equal at 1.8 x 10^ GV, and this is 
within less than a factor of 2 of the rigidities of both the proton peak and the iron peak. 
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and so the two propagation models will give almost identical fluxes for the same pulsar 
injection spectrum. The propagated fluxes are interestingly close to those observed, with 
the region between 3 x 10^ GeV and 2 x 10'^ GeV being about 20% of the observed flux. 

7 Discussion and Conclusion 

As can be seen in Fig. ^, we can expect fine structure in the spectrum in the region of 
the knee, as well as subtle composition changes. Our calculations were made assuming 
100% efficiency of conversion of spin-down power into accelerated iron nuclei, i.e. ^ = 1 
in Eq. |10|. Other uncertainties entering our calculation relate the propagation of cosmic 
rays. The diffusion coefficients assumed were based on an analysis of lower energy cosmic 
ray data, and their extrapolation to higher rigidities may be questionable. This would 
increase our predicted flux if the actual diffusion coefficient were lower. Also, the pulsar 
formation rate is rather poorly known. With all these uncertainties, we feel justified is 
adopting ^ ~ 1 in order see whether injection of cosmic rays accelerated in pulsar outer 
gaps could be potentially interesting. 

We shall next examine the composition in the knee region of the spectrum. The mass 
composition below ~ 10^^ eV is characterized by an average value of < In A >^ 1.2 ||^. 
If we assume that this component of the spectrum extends to higher energies, and add 
to it the contribution of heavy nuclei (with average A ^ 45) as obtained for our pulsar 
acceleration model, this would change the average value of < InA > to ~ 2.4 above the 
knee. Such a change in < In A > has recently been reported for the CASA-MIA air shower 
array data [Q. 

Our calculations have assumed a steady-state production of pulsar-accelerated cosmic 
rays throughout the Galaxy with a smooth spatial distribution of sources following that 
inferred for pulsars. It neglects the effect of relatively nearby sources, such as any within 
1 kpc releasing cosmic rays during the past 10^ years which would enhance the pulsar 
contribution, and so the effects of pulsar acceleration on the cosmic ray composition 
could be even more dramatic. Finally, we conclude that pulsar-accelerated cosmic rays 
could make a significant contribution of heavy nuclei to the cosmic rays observed in the 
knee region, and that this would manifest itself as fine structure in the spectrum, and a 
heavy composition at 3 x 10^-2 x 10'' GeV. 
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11 12 13 14 

\og{B sin 1 / G) 

Figure 1: Correlation between the minimum values of ellipticity k of the pulsars with age 
< 10^ years and their surface magnetic fields. The observed 7-ray pulsars are marked by 
crosses and other young pulsars are marked by stars. The initial periods of all pulsars 
have been fixed at Pq = 10 ms. 
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Figure 2: Dependence of the potential drop in the outer gap as a function of pulsar 
period for three values of surface magnetic field B = 10^^ G (dashed curve), 10^^-^ G 
(full curve), and 10^^ G (dotted curve). The corresponding thin lines show the maximum 
potential drop available across the region of the polar cap for three values of the magnetic 
field discussed above. For a given B value. Crab-type pulsars have the shortest periods. 
Vela-type pulsars have intermediate periods, and Geminga-type pulsars have the longest 
periods. 
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Figure 3: Radius of pulsar nebula (full line) as a function of age for a the pulsar with the 
parameters expected for the Crab pulsar and Nebula: surface magnetic field B = 3x 10^^, 
ellipticity A; ~ 3 x 10"^ 
supernova MgN = SMq, 



2^ , initial expansion velocity v = 2000 km s ^, mass of the 
and density of the surounding medium n = 0.3 cm~'^. The 



dashed line shows the radius, Rb, above which the magnetic field drops to the value 
of the interstellar magnetic field (assumed 5 x 10~^ G), and the dotted line shows the 
radius of the inner shock, i?sh, produced by the interaction of the pulsar wind with the 
surounding nebula. 
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Figure 4: Differential spectrum of nuclei injected by pulsars into the interstellar medium, 
multiplied by the square of energy per nucleus, (a) Contribution of different groups of 
nuclei: solid curve - A=l; dotted curve - A=2-10; short dashed - A=ll-20; dot-dashed 
- A=21-40; dot-dot-dot-dashed - A=41-56 (the contribution of A=56 is also shown 
separately as the long dashed curve), (b) Contribution of Crab (dashed curve), Vela 
(dot-dashed curve), and Geminga (dotted curve) type pulsars. 




Figure 5: Differential spectrum of nuclei injected by gamma-ray pulsars after propagation 
as described in the text (^ = 1). (a) D oc p^^^; (b) D oc p°-^^. Key to curves: thick solid 
curve - total; thin solid curve - A=l; dotted curve - A=2-10; short dashed - A= 11-20; 
dot-dashed - 21-40; dot-dot-dot-dashed - 41-56. The observed all-particle cosmic ray 
spectrum spectrum taken from ref. ^5 
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